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ABSTR^ACT 


This  thesis  is  a  study  of  a  computer  controlled  hydraulically  actuated  piston.  The 
syctem  uses  a  Hewlett  Packard  HP85B  microcomputer  as  a  controller.  Included  in  this 
research  is  a  detailed  computer  simulation  of  the  system  with  laboratory  validation. 
This  elTort  supports  the  overall  goal  of  complete  study  of  microcomputer  control  of 
electrohydraulic  power  systems  by  the  establishment  and  simulation  modeling  of  a 
baseline  system.  Special  emphasis  is  placed  on  modeling  the  effects  of  the  computer  on 
overall  system  performance.  It  was  found  that  sample  period  is  one  of  the  most 
important  factors  influencing  the  ability  to  control  a  hydraulic  power  element  using  a 
microcomputer.  Proper  selection  of  the  sampling  period  alone  is  not  always  sufficient 
to  insure  the  ability  to  control  the  plant.  Ocher  factors  such  as  non-linearities  in  the 
plant  may  influence  the  abilitv  to  use  a  digital  controller. 
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The  reader  is  cautioned  that  computer  programs  developed  in  this  research  may 
not  have  been  exercised  for  all  cases  of  interest.  While  every  effort  has  been  made, 
within  the  time  available,  to  ensure  that  the  programs  are  free  of  computational  and 
logic  errors,  they  cannot  be  considered  validated.  Any  application  of  these  programs 
without  additional  verification  is  at  the  risk  of  the  user. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Hydraulic  systems  in  one  form  or  another  have  been  around  for  a  long  time. 
These  systems  have  historically  relied  upon  humans  in  the  open  loop  case  or  analog 
electrical  circuits  in  the  closed  loop  case  for  system  control.  In  recent  years  however  a 
quiet  revolution  has  been  taking  place  in  the  hydraulic  systems  industry.  With  the 
current  rapid  advances  in  digital  technology  hydraulic  systems  are  fast  becoming  fertile 
ground  for  the  application  of  digital  control.  Not  only  are  digital  circuits  being 
designed  for  standard  control  functions  but  with  the  advent  of  microprocessors  these 
circuits  are  becoming  "thinking"  controllers.  It  is  this  trend  which  motivates  the 
detailed  study  of  microprocessor  controlled  hydraulic  systems. 

The  Department  of  .Mechanical  Engineering  has  begun  a  continuing  program  of 
study  on  the  subject  of  microprocessor  control  of  electrohydraulic  servomechanisms. 
The  overall  objective  of  the  work  is  to  provide  some  of  the  insights  necessary  for  the 
establishment  of  a  comprehensive  and  systematic  framework  embodying  the  design  of 
nucroprocessor-controlled  electrohydraulic  servomechanisms.  This  thesis  study  follows 
the  initial  work  done  by  Finch  [1]  in  an  effort  to  achieve  this  objective. 

1.  Advantages/Disadvantages  of  Hydraulic  Systems 

Hydraulic  systems  offer  several  unique  advantages  over  other  power  devices. 
These  systems  are  well  suited  for  applications  where  low  component  weight,  high 
accuracy,  and  heasw  load  movement  is  required.  .Merritt  (2)  states  the  following 
specific  advantages  of  hydraulic  systems; 

1.  Smaller,  lighter  more  compact  components. 

2.  Hydraulic  fluid  acts  as  a  lubricant  which  e.xtends  component  life. 

3.  There  is  no  phenomenon  in  hydraulic  components  comparable  to  the  saturation 
and  losses  in  magnetic  materials  of  electrical  machines. 

4.  Hydraulic  actuators  have  a  higher  speed  of  response  with  fast  starts,  stops,  and 
speed  reversals  possible. 

5.  llvdraulic  actuators  may  be  operated  under  continuous,  intermittent,  reiersing. 

■  ir.d  'tailed  conditions  without  damage. 

6.  Hydraulic  actuators  have  a  higher  stillness  so  there  is  little  drop  in  speed  as 
loads  are  applied. 
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7.  Open  and  closed  loop  control  of  hydraulic  actuators  is  relatively  simple  using 
valves  and  pumps. 

There  are  some  disadvantages  to  the  use  of  hydraulic  systems  and  no 
discussion  would  be  complete  without  at  least  mentioning  them.  .Merritt  [2j  lists  these 
shortcomings  as  follows: 

1.  Hydraulic  power  is  not  as  conveniently  available  as  electrical  power. 

2.  Small  allowable  tolerances  result  in  high  costs  for  the  manufacture  of  hydraulic 
components. 

3.  The  hydraulic  fluid  imposes  an  upper  temperature  limit. 

4.  Hydraulic  fluid  may  present  a  fire  hazard. 

5.  Hydraulic  systems  are  messy. 

6.  Hydraulic  systems  are  susceptible  to  failure  due  to  dirt  or  contamination  of  the 
fiuid. 

Although  these  disadvantages  may  preclude  the  use  of  a  hydraulic  system  in 
some  specific  applications,  their  widespread  use  attests  to  their  versatility  and 
usefulness.  Hydraulic  systems  are  of  the  most  use  in  applications  where  relatively  high 
power  levels  are  required  especially  within  confined  spaces.  This  explains  their 
widespread  use  in  aircraft  control  surface  positioning,  fin  stabilizer  positioning  on 
ships,  and  even  as  a  means  of  propulsion  in  special  applications  such  as  on  some 
SW'ATH  vessels. 

2.  .\dvantages  of  Digital  Microprocessor  Control 

Until  recently  hydraulic  systems  have  relied  upon  analog  circuitry  for  control. 
.More  and  more  these  systems  are  being  manufactured  with  digital  control.  One  would 
logically  ask  what  advantages  this  offers  over  the  proven  analog  controller.  Aside  from 
the  increase  in  reliability  and  decrease  in  signal  voltages  necessary,  the  major  advantage 
to  digital  control  is  the  ability  to  incorporate  a  microprocessor  in  the  circuit.  This 
improvement  broadens  the  range  of  functions  that  may  be  included  in  the  control  loop. 
Henke  [3|  lists  the  following  new  functions  for  the  controller: 

1.  Closing  the  loop,  replacing  the  analog  summing  point  and  conditioning  the 
error  signal. 

2.  Pre-loop  processing 

3.  Peripheral  Processing 

4.  .Adaptive  Control 

5.  Smart  Redundancy 

6.  Time  Optimal  Control 
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Along  with  increasing  the  range  of  applications  there  are  other  major 
advantages  to  the  use  of  digital  microprocessor  control  circuits  for  hydraulic  systems. 
No  doubt  a  leading  advantage  is  cost.  As  the  complexity  of  the  control  system  rises 
analog  control  circuitrx-  costs  rise  dramatically.  Although  this  trend  is  seen  with  digital 
systems  the  rise  in  cost  is  less  dramatic.  In  addition,  the  cost  of  changing  an  analog 
control  circuit  once  in  place  is  much  greater  than  a  digital  circuit  which  incorporates  a 
microprocessor  [4]. 

Another  major  advantage  of  digital  microprocessor  control  is  flexibilitv. 
While  analog  circuits  may  require  component  changes  and  rewiring,  digital 
niicroprocessor  circuits  may  require  only  reprogramming.  This  may  be  accomplished 
by  either  keyboard  entry  or  by  simply  reprogramming  a  memorx'  chip.  This  reduces 
the  time  and  cost  of  making  changes  to  the  control  system.  It  may  also  enable  the 
overall  system  to  be  used  for  several  functions  rather  then  one  specialized  task. 

3.  Previous  Thesis  Work 

In  Ref  1.  Finch  reports  the  results  of  his  study  of  a  microcomputer  controlled 
hydraulic  motor.  Me  used  a  512K  Macintosh  computer  and  a  product  of  GW 
Instruntents  called  Mac.ADlOS  as  as  interface  to  control  a  hydraulic  system  consisting 
of  a  hydraulic  power  source,  servovalve  and  hydraulic  motor.  In  this  work  the 
computer  was  used  for  both  open  and  closed  loop  signal  processing  and  conditioning. 
Finch  was  successful  in  implementing  the  computer  as  a  system  controller  and 
establi-hing  a  baseline  sy.stem  for  follow-on  work.  He  used  Dynamic  Simulation 
Language  (DSLi  programs  to  predict  system  performance  but  his  analysis  did  not 
include  any  digital  microprocessor  dynamics  in  these  simulations. 

B.  OBJECTIVES 

Finch,  made  several  recommendations  for  follow  on  work.  One  of  these  wa';  that 
piston  actuator  systems  be  studied  in  further  detail  since  "a  large  varietv  of  hvdraulic 
‘systems  involve  the  positioning  of  loads  using  a  piston-type  arrangement."  One  of  the 
main  objectives  of  this  thesis  cfTort  has  been  to  study  the  hydraulically  controlled 
piston  response  in  detail.  .A  Hewlett  Packard  HPS5B  microcomputer  has  been  used  as 
the  system  controller.  .An  additional  objective  was  to  model  the  computer  dvnamics 
for  a  more  accurate  sim.ulation  of  the  overall  .svstein. 
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C.  APPROACH 


To  begin  the  study  a  thorough  understanding  of  the  analog  system  was  sought. 
This  began  with  a  disection  of  the  system  into  individual  components  so  that  a  linear 
mathematical  model  could  be  developed  for  each.  The  component  equations  were  then 
cembined  and  coded  in  a  simulation  using  DSL.  Step  input  performance  predicted  by 
the  simulation  was  verified  by  lab  testing.  Once  the  analog  system  was  fully 
understood  and  confidence  was  gained  with  the  simulation  a  microcomputer  control 
system  was  added.  Adding  this  system  required  the  construction  of  a  digital  or  discrete 
time  model  to  represent  it.  This  was  then  combined  with  the  analog  model  to  develop 
the  overall  disital  svstem  simulation. 


coils.  The  power  amplifier  gain  may  be  varied  from  110  to  1520  nrLA  volt.  Figure  (2.3) 
shows  the  relationship  between  current  output  and  voltage  input  of  the  servoamplifier 
at  the  two  extremes  of  the  amplifier  gain. 


B.  COMPONENT  EQUATIONS 
1.  Seno Valve  Dviiamics 

The  servovalve  used  was  found  through  experimentation  to  be  slightly  open 
centered  but  will  be  considered  critically  centered  for  the  purposes  of  modeling.  (The 
degree  of  open  centered  behavior  is  so  slight  that  the  more  complicated  model  is  not 
warranted.)  The  general  equation  describing  the  operation  of  a  critically  centered 
spool  valve  is: 


flow  through  the  load 
discharge  coefficient 
area  gradient 
valve  spool  displacement 
mass  density  of  the  hydraulic  fluid 
supply  pressure 
load  pressure  drop 


.Modification  to  Eq.  (2.1)  is  necessary  due  to  the  difficulty  in  directly 
measuring  valve  spool  movement.  To  reflect  this,  the  following  parameter  which 
relates  valve  spool  movement  to  current  applied  to  the  valve  torque-motor  coils,  is 
defined; 


C  ,w 
G 


^/z 


(P) 


(2.2) 


With  this  substitution,  Eq.  (2.1)  becomes: 


(:.3) 
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K^.  =  valve  current  transfer  coefficient 
(3.52xlO'^:n^  sec-ma-psi^  *■) 

I  =  current  supplied  to  the  valve 

Note  that  the  use  of  the  modified  valve  current  transfer  coefficient  presumes  a  linear 
relationship  between  the  current  supplied  to  the  valve  coils  and  valve  spool  position. 
Experiments  have  proven  this  to  be  a  valid  assumption. 

Equation  (2.3)  must  now  be  linearized  to  facilitate  modeling  and  studies  of 
system  response.  The  no-load  null  valve  spool  position  is  chosen  as  the  center  point 
for  the  lineanzation.  Reasons  for  this  choice  are  "system  operation  usually  occurs  near 
this  region,  the  valve  flow  gain  is  largest,  giving  high  system  gain,  and  the  flow- 
pressure  coefficient  is  smallest,  giving  a  low  damping  ratio.  Hence  this  operating  point 
is  the  most  cntical  from  a  stability  viewpoint,  and  a  system  stable  at  this  point  is 
usually  stable  at  all  operation  points.''[21  The  linearized  form  of  Eq.  (  2.3)  is: 

QL=Kqi^^-^ci^^L  (2.4) 


where 


K  .=K  (P^-P. 
qi  V  s  L 


Kq-  =  modified  null  flow  gain 


0.5  K^AI 


K  =  modified  null  flow  rressure  coefficient 


Subscr.pt  0  refers  to  conditions  at  the  operating  point  and  a  4  denotes  departures  from 
this  point.  Since  we  are  linearizing  about  the  null  position,  i.e.  Q^q  =0,  Iq  =u,  and 
Pr  Q  =0,  Eq.  (2.4)  may  be  written: 
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2.  Piston  Dynamics 


For  the  piston  actuator  the  model  suggested  by  Merritt  [2]  is  adopted  as 


follows: 


Q  =A  X  +C^  P_  + 
p  p  tp  L 


where  the  (•)  notation  indicates  rate  of  change  and 
Ap  =  piston  area 

Cjp  =  total  leakage  coefficient  of  piston 

=  total  volume  of  piston  chamber  and  supply  hose 
e  elfective  fluid  bulk  modulus 


The  total  leakaae  coefficient  is  given  bv: 


where 


C  =C.  +0.5C^^ 
tp  iP  ep 


C-p  =  internal  leakage  coefficient 
=  external  leakaae  coefficient 


For  the  system  under  investigation  in  this  study,  the  external  leakage  has  been 
deterrruned  to  be  negligible  and  Cjp  =  C-p.  With  this  simplification,  and  after  Laplace 
transformation  (with  zero  initial  conditions),  Eq.  (2.6)  becomes: 


(s)=A  sx^(s)+C.^P.  (s)+C^^sPt  (s) 
L  P  P  tp  ii  Cp  L 


where  the  I'luid  caoacitance  is 


C  =  T 

cp  4 


3.  Load  Dynamics 


The  load  is  modeled  as  a  simple  spring,  mass,  dam.per  system.  In  the  general 


F  =M  X  +3  Ik  +KX  +F. 
g  t  p  P  P  P  L 


=  force  developed  by  piston 
=  total  mass  accelerated 
=  VISCOUS  load  damping  coefficient 
=  spring  load  constant 
=  arbitrarv'  load  force  on  piston 


The  force  generated  by  the  piston  is  determined  by  the  piston  area,  the  load 
pressure  drop,  and  an  elTiciency  to  account  for  losses  in  the  actuator.  In  this  specific 
application  there  were  no  spring  forces  so  K  =  0.  The  equation  is  now  transformed  into 
the  Larlace  domain  with  the  the  following  result: 


( s )  =M^s  ^  ( s )  +BpSXp  ( s ) 


C.  ANALOG  SYSTEM  MODEL 

Now  that  the  linear  mathematical  model  for  each  component  of  the  system  has 
been  developed,  these  equations  can  be  combined  to  formulate  the  model  for  the 
hydraulic  power  elem.ent  with  servovalve  current  as  an  input.  First  Eq.  (2.5)  is  equated 
to  Eq.  (2.1')  to  eliminate  Qj_.  This  combined  equation  is  then  rearranged  to  express 
as  a  function  of  the  input  current  and  the  piston  velocity.  The  result  is: 

(K  .  I-A_sx_  (s)  )  (= - — )  /Tin, 


K  . +C.  ^ 

Cl  ip 


rearrangement.  Eu.  (2.9)  m.av  be  written: 


.  a  D  { s )  )  (— — ) 


sx  (3)=:<  (s)=- 

Z)  o 


3  .  1) 

P 


Equations  (2.10  i  and  i  2.1 1 )  >:u2gcst  the  defnition  of  two  time  constants,  as  follows 


*el 


*'.<3 


SI 

% 


>sa 


\%here  ^ci'^^ip  leakage  coefficient).  The  resulting  block  diagram  is 

shown  in  Ficure  <  2.4). 


Figure  2.4  .Analog  System  Block  Diagram. 

D.  analysis  of  the  system 

From  preliminary  lab  observations  the  closed-loop  analog  svstem  exhibits  Hrst 
order  response  behavior  at  low  frequencies.  It  is  clear  from  Fig.  (2.4)  that  the 
hsdrauiic  power  element  is  third  order.  In  an  effort  to  understand  this  apparent 
..ontradiction  lurther  analytic  analysis  of  the  system  is  necessarv.  The  results  will  <er.  e 
as  a  ber.cnmark  for  the  computer  simulation  and  system  expcnmentation.  Of  prime 
nterc-t  are  the  system  time  constants  and  break  frequencies. 

Fi.-st  some  block  diagram  reduction  is  useful.  Under  no-load  '  t:.e  mst.-n 

■  elcci'y  may  be  represented  by  the  following  transfer  function: 

sx  (s)  -A 


lab 


With  some  further  simplification  Eq.  (2.12)  may  be  rewritten: 


G  (.s).= 
P 


(2.13) 


where: 


B  K 
o  ce 
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p  ce 
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Ti  +  T2 


B  K^' 
o  ce 
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Recall  incorporates  the  piston  internal  leakage  (Cjp)  and  the  modified  null 
few  pressure  ^cetEcient  (Kj^j)-  It  is  known  from  e.xpenmentation  and  observation  of 
the  s'.stem  that  Cjp  is  on  the  order  of  10'*'  in  sec-psi.  What  is  necessaip'  is  to 
e'tin'a'c  the  range  of  values  over  which  will  var>'.  Referring  to  Eu.  i  the 
rr.  !\.r;t^rr.  value  of  K  ,,  will  occur  when  I  is  Z'tt'  rrL-X  and  the  valve  ;s  ;'u!;\  oren. 
Pu’t.ng  ;r,  the  ’.aiues  for  this  appucanen  iPj,  =  S<Mj  psi.  =  obtain  a  value  of 

'1012  in  '  sec-psi.  A  lower  bound  must  now  be  found,  .\ccording  to  the  servovai'.e 
manufacturer  the  maximum  quiescent  flow  at  the  null  valve  position  is  0.3.5  gpm  at 
'"'ll)  psi.  The  system  used  for  experimentation  was  run  at  SiW  psi  so  one  would  expect 
tiiat  U'lng  the  jOoo  psi  quiescent  How  figure  would  give  a  conservative  estimate  of  the 
hr-ver  ‘•'o.ind.  L  M.ng  thus  figure  for  the  lower  bound  the  overall  range  of  values  for  K^_.. 


iA 


mim 


is  from  5.0x10"^  in^;lb-psi  to  1.2x10'^  in^;lb-psi.  In  response  to  step  inputs,  the  valve 
is  fully  open  for  a  majority  of  the  time  so  that  a  value  of  l.4xl0''^  in-’  Ib-psi  has  been 
selected  for  under  these  conditions. 

Noting  that  K^g  and  Bp  are  very  small  so  that  division  by  their  product  is  a 
number  which  is  very  much  greater  than  I,  the  following  simplifications  to  the 


constants  of  Eq.  (2.13)  are  made: 


ce 


The  final  block  diagram  of  the  closed-loop  system  is  shown  in  Fig.  (2.5),  with  Ga 
denoting  the  power  amplifier  transfer  function  (assumed  to  be  a  simple  gain  in  this 
analysis). 


Figure  2.5  Closed  Loop  .Analog  .System  Block  Diagram. 
From  this  block  diagram  the  overall  t.'-ansfcr  function  is  given  as: 


When  doing  laboraton-  experiments  the  output  (Xp)  is  compared  to  the  input  (e^) 
in  its  teedback  form  Cp.  Noting  that  ep(s)=  H(s)Xp(s)  and  making  the  substitution  for 
the  plant  transfer  function  the  closed-loop  transfer  function  becomes: 


e_  (s) 


^  s  -t-  s  +  1)  +  K 


(2.15) 


!P 


I 


where 


K=  G  (s)K  H(s)K  . 
A  p  qi 


In  order  to  make  an  estimate  of  the  break  frequencies  of  the  closed  loop  system 
some  work  with  the  characteristic  equation  is  required.  The  characteristic  function  of 
Eq.  1,2.15)  may  be  approximated  as: 


under  the  condition  that 


(s+K)  [s  +(2Cj^cJj^-K)  S  +  Uh  ] 


2Ch%-!'='0 


This  condition  is  satisfied  in  the  present  case  and  Eq.  (2.15)  may  therefore  be 
approximated  as: 


e  (s) 

p _ 

e^TsT 


!s-t-K)  (s'+'Ji 


(2.1^' 


The  values  for  K  and  arc  ?0  sec'^  and  261  sec'^  respectively  i^see  Appendix  A  I'or 
constant  values  used  to  obtain  these  results).  This  yields  a  break  frequency  and 
hydraulic  natural  frequency  of  about  4.S  Hz  and  dl.6  Hz.  respectively. 


!■ 

IThe  tonclucicn  to  be  drawn  from  this  analysis  is  that  the  analog  closed  loop 

system  acts  like  a  first  order  system  for  low  frequencies  and  a  third  order  system  with 
nunimal  damping  at  the  higher  frequencies.  The  hydraulic  power  element  has  very  low 
damping  which  would  make  it  unstable  were  it  not  for  the  integrating  effects  of  the 
'.s- K't  term  in  the  characteristic  equation. 


E.  FREQUENCY  RESPONSE  TEST 

To  verify  the  analytic  results  of  the  previous  section,  frequency  response  tests 
were  conducted.  The  results  are  shown  in  Fig.  >2.6a  and  2.6b).  From  these  figures  two 
break  frequencies  are  observed  within  the  rar.ge  of  measurable  system  output,  one  at 
2."  Hz  and  one  at  ITS  H/.  The  second  break  frequency  lies  near  the  value  estimated 
for  the  ■'Crvovalve.  which  heretofore  has  been  modeled  as  a  simple  gain. 

.Although  the  frst  break  frequency  of  2.”  Hz.  is  less  than  the  predicted  value  of 
-l.S.  this  can  be  explained,  at-least  partially,  by  the  uncertainty  in  the  values  leading  to 
this  estimate.  In  addition,  the  amplitude  of  the  input  sine  wave  in  these  tests  caused 
the  v.tl\c  to  be  actuated  over  43  of  its  range  of  motion  to  either  side  of  null.  .Such 
l.trgc  excursions  may  be  outside  of  the  region  where  the  valve  flow  linearization  is 
valid.  With  the  result  that  the  amplitude  plot  is  attenuated  below  the  predicted  linear 
response.  In  any  case,  the  predicted  dominance  of  the  first  break  frequency  is  verified 
by  the  tests, 

F.  MODEL  VALID.ATION 

Once  a  thorough  theoretical  understanding  of  the  analog  system  was  achieved  a 
DSL  simulation  was  written  using  Fig.  (2.4)  and  Fig.  {2.>)  as  guides,  fho  simulation 
code  is  given  ;n  .Appendix  .A.  Figure  (2.'}  shows  a  comparison  of  experimentally 
obtained  data  and  the  model  ,^or  a  0.1^5  in.  step  input.  .Agreement  between  predicted 
and  actual  performance  is  good  at  low  step  inputs  however  as  the  size  of  the  step 
increases  the  model  prediction  becomes  less  accurate.  Figure  (2.S)  shows  a  l.'.S.s  in. 
step  input  response.  Step  response  comparisons  for  steps  between  these  values  are 
given  in  .Appendix  B. 

The  most  likely  explanation  for  the  degradation  of  model  performance  as  tb.e  step 
size  increases  has  to  do  with  the  valve  linearization.  Since  this  linearization  was  done 
around  the  null  valve  position  one  would  expect  more  model  prediction  error  with 
increasing  step  ^ize.  This  is  because  as  the  step  size  is  increased  large  excursions  of  t'ne 
valve  .ire  reciuired.  Fventually  the  step  size  becomes  so  large  that  the  valve  goes  to  the 
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Figure  2.6a  Closed  Loop  Frequency  Response  (Amplitude) 
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Figure  2.ob  Closed  Loop  Frequency  Response  (Phase) 


fully  open  position.  Once  this  occurs  increasing  the  step  size  further  causes  the  valve 
to  remain  in  this  fully  open  position  longer. 

Still  another  explanation  lies  in  the  values  of  some  of  the  constants  used.  The 
model  is  sensitive  to  changes  in  the  load  damping  constant  (Bpj  which  was  obtained 
empirically.  The  model  also  shows  sensitivity  to  the  value  of  the  modified  null  flow 
coelTicient  (K^;).  Both  of  these  constants  were  refined  from  initial  estimates  to  obtain 
the  model  performance  shown  in  Fig.  i:.?!  and  Fig.  (2.8).  Similar  performance  may 
possibly  be  obtained  by  the  variation  of  the  actuator  efficiency  and  the  value  for 
the  effective  bulk  modulus  (S^  which  are  also  empirically  based. 


III.  DIGIT.4L  SYSTEM 


A.  DESCRIPTION  OF  THE  SYSTEM 

The  digital  system  includes  all  the  components  of  the  analog  system  plus  the 
following  Hewlett-Packard  computer  and  interface  equipment  which  together  comprise 
the  microcom.puter  control  system. 

HP  85B  computer 

HP  S2901M  twin  flexible  disc  drive 

HP  7470A  plotter 

HP  6942A  multiprogrammer  with: 

2  A  D  converter  cards 
1  D  A  converter  card 
1  high  speed  memory'  card 
1  memory  expansion  card 
1  timer  pacer  card 

Figure  (3.1)  shows  the  complete  digital  system.  Where  previously  the  reference 
signal  was  fed  to  the  power  stage  input  (terminal  4)  on  the  servoamplifier  the  error 
signal  from  the  digital  controller  is  now  applied.  The  feedback  signal  is  fed  to  the 
digital  controller  so  that  the  dilTerencing  operation  is  shifted  from  the  servoamplifier  to 
the  microcomputer  controller. 

The  operation  of  the  system  under  study  begins  with  the  operator  entering  a 
desired  position  in  inches  at  the  microcomputer  keyboard.  At  this  point  the  computer 
algorithm  (see  .Appendix  C)  samples  the  feedback  signal  to  determine  the  load  position. 
This  information  is  passed  through  the  controller  portion  of  the  algorithm  and  an  error 
signal  is  generated  and  sent  to  the  servoamplifer.  The  generation  of  the  error  signal 
continues  until  the  load  reaches  the  desired  position.  The  algorithm  includes  a  position 
cliecking  routine  which  maintains  the  commanded  position  until  a  new  value  is  entered 
at  the  keyboard. 


of  plant  control  used.  The  sampling  rate  is  not  merely  the  time  the  analog-to-digital 
(A;D)  and  digital-to-analog  (D'A)  converters  take  to  process  signals.  It  must  also 
include  the  time  taken  for  the  controller  algorithm  to  process  the  input  signals  and 
generate  an  output  signal.  Thus  the  sampling  period  "T"  is  the  total  time  which 
elapses  between  sampling  of  the  feedback  signal  and  transmittal  of  the  error  signal. 

The  elTects  of  sampling  are  more  clearly  understood  if  first  an  example  of  how 
the  system  operates  is  presented.  Let  us  say  the  load  is  at  the  center  (null)  position. 
.At  time  t  =  0  a  1-in.  position  is  requested  at  the  keyboard.  This  causes  the  position 
feedback  input  to  be  sampled,  inverted  and  summed  with  the  requested  position.  Next 
the  summed  signal  is  passed  through  the  controller  algorithm  and  the  resulting  error 
signal  transmitted  to  the  servoamplifier.  The  error  signal  remains  at  this  level  until  one 
sampling  period  (IT)  has  elapsed  and  it  is  updated.  For  example,  if  the  sampling 
period  were  100  ms  and  the  inital  error  signal  10  volts,  10  volts  would  be  continuously 
sent  to  the  serx'oamplifier  for  a  period  of  100  ms.  Then  presumably  the  error  signal 
would  be  updated  to  say  5  volts  and  5  volts  would  be  continuously  sent  to  the 
servoamplifier  for  the  next  100  ms  etc. 

In  the  analog  system  a  smooth  continuous  error  signal  resulted  from  the 
summation  of  the  reference  and  position  feedback  signals.  In  the  microcomputer 
control  system  the  error  signal  produced  is  a  continuous  series  of  steps.  .A  signal  of 
this  form  is  said  to  be  "quantized."  The  steps  have  a  duration  of  one  sample  period.  If 
the  sample  period  is  too  long  then  the  error  signal  is  not  changing  fast  enough  to 
control  the  system.  This  will  cause  the  system  to  be  unstable.  To  prevent  this,  a  rule 
of  thumb  adapted  from  Shannon's  Sampling  Theorem  is  used.  The  rule  suggests  that 
the  sampling  frequency  be  ten  times  as  large  as  the  highest  break  frequency  of  the 
system. 

One  of  the  ways  to  model  the  behavior  of  the  microcomputer  control  system 
is  shewn  in  Fig.  (3.2).  While  there  are  several  types  of  holding  devices  the  zero-order 
hold  ;s  best  suited  for  this  application.  The  zero-order  hold  takes  in  a  digital  value  and 
outputs  that  value  continuously  until  one  time  period  has  passed  and  it  receives  an 
updated  digital  value.  The  input  to  the  zero-order  hold  is  e^(t)  Jwherc  the  (")  denotes  a 
digitai  or  discrete  time  signal}  which  is  the  error  signal  generated  by  passing  the 
summed  signal  e^(t)  through  the  controller  C^.  The  output  from  the  zero-order  hold 
is  e/kT)  the  quantized  error  signal.  The  combination  of  the  switch  and  the  zero-order 
hold  represent  the  D  .A  converter. 


Figure  3.2  Microcomputer  Controller. 


The  plant  controller  represents  that  part  of  the  algorithm  which  contains 
the  controller  logic.  This  logic  can  take  the  form  of  a  simple  gain  or  be  more 
complicated  such  as  a  proportional-integral-derivative  controller.  To  model  this  a 
discrete  time  representation  of  must  obtained.  To  help  illustrate  the  process  an 
example  is  presented.  The  approach  begins  by  obtaining  the  Laplace  domain 
representation  of  the  controller  transfer  function  {G^(s)}.  Say  the  controller  uses 
integral  control  action.  This  is  represented  in  the  Laplace  domain  as: 

e*(s)  K.  (3.1) 

e  _ 


Gc(s)  - 


Next  G^(s)  must  be  transformed  to  the  Z  or  discrete  time  domain  by  obtaining  the  Z- 
transform.  This  may  be  done  by  one  of  three  methods:  by  transforming  first  to  the 
time  domain  and  then  to  the  Z-domain  (the  rigorous  approach),  by  an  approximation 
method  such  as  Euler's  or  backward  dilTercnce.  or  by  using  tables.  For  this  example  a 
table  was  used  to  obtain  the  following: 


G(z) 


e  (z) 
e 


Dividing  the  top  and  bottom  of  Eq.  (3.2)  by  z  the  following  results: 


z 


With  some  rearrangement  Eq.  (3.3)  can  be  wntten: 


e  (z)  =  K.e  (z)  +  e  (z)z”^ 
e  1  e  e 


The  inverse  Z-transform  is  now  taken  for  each  term  of  Eq.  (3.4)  to  obtain  the  discrete 


time  mouei. 


e*  (t)  =  K^e*'  (t)  +  e*(t-T) 


The  e4(t-T)  term  represents  the  value  of  ea(t)  one  time  period  backwards  and  is  knowT 
as  the  lirst  past  value  of  e_,(t).  Had  the  e.(z)  been  multiplied  by  z  in  Eq.  (3.5)  then 
this  would  have  inversely  transformed  to  e.it-nT)  which  represents  the  value  of  e^(t)  n 
time  penods  in  the  past. 

Had  Eq.  (3.3)  been  more  complex  a  more  ngorous  approach  would  have  been 
employed.  The  equation  vvould  have  been  arranged  in  such  a  manner  as  to  permit  its 
representation  as  a  geometnc  series  or  a  combination  of  several  geometric  senes.  Then 
the  remainder  of  the  procedure  would  be  followed  as  in  the  previous  example  with  the 
complication  that  rather  than  a  finite  number  of  terms  the  number  of  terms  I'in  Eq 
3.5)  far  exam.ple!  would  be  infinite.  A  decision  would  have  to  be  miade  as  to  how 
many  terms  to  use  to  approximate  the  series.  Usually  only  the  first  few  terms  are 
needed  to  yield  suiTcient  accuracy  for  modeling. ['.(,') 

For  this  specific  application  a  proportional  controller  is  used  s-  -ii,;; 
In  this  case  it  is  not  necessan.'  to  •lo  throuch  the  various  transformati-sns 

e  pi.  .  ..  _ 

Since  K..., ,  is  merely  a  constant  and  remains  so  in  the  discrete  time  Jom.ain. 

pc 

C.  DIGIT.\L  SVSTE.M  MODEL 

The  digital  system  m.odel  is  a  combination  of  the  continuous  time  plant  model 
deveicned  in  chapter  two  and  the  discrete  time  rrucrocomputcr  controller  model 


explained  in  this  chapter.  The  block  diagram  of  the  system  model  is  shoum  in  Fig, 

^  ^  \ 

I  ). 


Figure  3.3  Digital  Svstem  Model. 


D.  MODEL  V.ALIDATION 

Once  conl'idence  was  gained  in  the  validated  analog  simulation,  coding  was  added 
to  model  the  microcomputer  control  system.  The  biggest  problem  in  obtaining  an 
accurate  simulation  was  predicting  the  value  for  T.  A  value  of  0.35  sec  was  chosen 
based  on  timing  the  execution  of  the  controller  segment  of  the  computer  algorithm. 
Using  this  value  and  Kp^=  1.0  the  simulation  predicted  a  stable  limit  cycle  response  to 
any  size  step  input.  With  K,.,j,  =  0.05  the  simulation  predicted  stable  responses  to 
■.aricus  step  inputs  simalar  to  those  obtained  from  the  analog  system.  .At  this  point  lab 
validaticn  of  these  results  was  attempted. 

The  prediction  of  performance  with  =  !.<■)  was  qualitativelv  correct.  Stable 
'..I'.'.it  cle  performance  evas  observed  regardless  of  the  size  of  the  ''tep  input.  Next 
K.,.=  >.'is  was  attempted.  .\t  first  :t  appeared  that  again  the  simulation  had  preuicted 
tn-e  perfcrmance  accurately.  ,-\s  time  wore  on  however  step  responses  became  mere 
and  .more  oscillators'  untii  once  again,  stable  limit  cycle  response  to  anv  size  step  was 
or'C’". ed,  .A;'rer  some  though:  it  became  clear  that  '.viiat  had  varied  o'.  cr  tb.e  tim.e  of 
tl'.e  testing  was  the  temperature  of  the  oil.  .-\t  the  beginning  ol'  the  tests  the  oil 
tentrerature  was  between  bo  F  and  ""f  1-  and  stable  step  response  beiiavior  was 
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observed.  .As  the  oi!  heated,  mare  and  more  oscilIator>’  response  was  obscr'. eoi 
Finady,  around  S5°F  limit  cycle  response  resulted.  With  these  findings  an  operat.n. 
emperatare  range  of  lOil'^F  to  !  lOT  was  chosen  for  the  remainder  of  lab  tcstinc.  1;..' 
s  not  to  say  that  the  failure  of  the  microcomputer  to  control  the  plant  .s  .cleiv  a  resul: 
if '.a.nations  in  the  oil  temperature.  In  fact  the  analog  system  was  stable  rega'dless  oi 
he  'ernperature  and  one  would  expect  similar  performance  from  the  digital  sysicm. 

The  system  performance  with  regard  to  oil  temperature  is  interesting  to  note,  B;. 
vay  of  explanatio.a  of  this  behavior  as  the  tem.perature  goes  up  the  oil  iscositv 
lecreases  by  a  factor  of  about  10  in  the  range  6b“F  to  100° F.  Tins  in  turn  decreases 
he  VISCOUS  damping  on  the  already  highly  underdamped  plant.  This  is  true  .''or  both 
he  analog  and  digital  systems,  however,  and  still  dees  not  explain  the  failure  of  me 
ligital  simu.ation  to  accurately  predict  system  performance,  given  an  aue.juate 
ampiing  I'requency. 

One  of  the  cr.tica!  parameters  for  stable  digital  control  is  the  sample  period. 
,  'ing  the  rule  of  thumb  suggested  by  Shannon  s  sampling  thecremi  a  sampling  period 
f  no  mere  than  Z.'fxli)'-  sec  should  be  used.  This  value  is  obtained  using  the  '.afue  oi 
he  hvdrau.ic  natural  frequency  as  the  highest  plant  break  frequency.  This  insures  an 
dequate  sampling  per.cd  but  does  not  account  for  e.Tects  of  changing  the  ..or.tro.ler 
am.  Reducing  the  controller  gam  slows  the  response  time  of  the  entire  s;. stem,  and 
hus  a  longer  sampling  period  may  be  tolerated.  To  refect  this  a  slight  variation  o!  the 
ule  of  thumb  is  used.  Rather  than  using  ten  times  the  highest  break  frequenc.',  ten 
imes  the  open  loop  bandwidth  of  the  system  frequency  response  is  empicyeu,  liie 
and'.vidth  IS  approximately  K.  the  open  loon  gain.  With  1  "  a  sample  period  of 

'xl'i''  sec  was  obta.ned.  This  means  our  controller  is  sampling  too  sloub.  to  ^.mtr'l 
ne  plant.  \\'itii  =  0.'''5  a  sample  period  of  0.5  sec  is  obtained.  Tiii';  mdlca•e^  •f.at 

oe  ..iigiml  .ontrolier  using  this  gam  sheuid  be  able  to  control  tite  plan’,  I  hi^  '-er.e^  •  ' 
-"'taiim.te  me  predictions  made  bv  tite  simulation  but  doec  not  iielp  exr.,;::i  .^.o  me 
'  'em  e■.;l;''l•^  limi'  cvci.e  performance  e\en  I'ltii  me  lower  ..)nT.'.ler  con 


.n'- e^tioateu  was 


. .  .  k.  . 


.v'L*: 


:is  :e. 


^ '  L' 


::;e  ■.  :!;ucs  :cr  :ne  poniiior;  tne  x;urrc:n  ^upri.ed  ro  the  serve., n’-e  1.,  :: 
Tor  Munal  '.alues  were  observed.  .A  n^'men  of  the  data  run  is  mc.uued  m  A-r 


I\  .  C  ONCLL  SIGNS  AND  RECOMMENDATIONS  FOR  FOLLOW  ON 

WORK 


.'V 


f, 

m 


K.  (  ()N(  LI  SIGNS 

\1  ro.  .':r.p  ..;cr  ^j:uro!  of  h\Jrai!!uaii\  aLtuateJ  steins  is  possible  ansi  oilers 
mans  aJ'- ..nt  iees  o\cr  s  am enti anal  analce  control.  Of  estreine  iir.portance  in  sligitai 
^  \  •.'■'.e  vari'.ple  peris'-J.  Selection  of  control  harduare  soltuare  siiouid  nsrt  be 


unacr'.ir.en  itn 


first  ui.slersrai.dine  tlie  plant.  .Xn  estimate  of  tl'.e  h.ei'.est  brea.k 


: . e  , i.e:.-'.  :;'...st  'e  obtained  to  determine  tiie  necessary  sampling  period.  OnK  alter 
t;.;-  .s  determined  can  one  properlc  select  hardware  software  for  the  plant  controller. 
Ht.t  proper  seiestion  i.f  tne  sampling  period  is  not  enough.  Non-Iinearities  m  tiie  plant 
mat  >.  case  .r.stabi.ities  at  the  dagitalis  controlled  ssstem.  It  is  theon/ed  by  tf.c 
n.sc  under  nh.!:  fnrtlicr  reduvtions  ai  tl.e  sample  period  will  osercome  this  problem  for 
•fe  S-,  xwm  staviicd.  It  .'  as  not  possible  to  serif’s  this,  hut  the  last  tltat  th.e  analog 
'  .••em  s'a's.e  -.ert .inti's  ntdwnitcs  tiiat  somewhere  between  .1  santple  period  of 
ssw  .;nd  ''''  'sn  statsOm.  js  i.nic'.ed. 


B.  RECOMMENDATIONS  FOR  FOLLOW  ON  WORK 
!  \  al\e  Linearization 

( )ne  nf  :f.e  greatest  '.seaknesses  0!' the  analog  and  digital  models  he  'svitii  ’he 
.  nc.ii.  n  :  the  saise  Idns  equa'-nit.  Miwh  ot  the  time  the  sad-e  w  opcr.nmg  ./  tne 
■  I'cn  p’'s;-;,’i:.  r.ither  titan  mad  position.  Because  a  choue  0!' one  posonen  iiad  '0 
w  1C  : .  '  'ttesis  rescarcit  rite  null  position  rsas  mes'  .ogicai  lor  reasons  airc.idr 
n  jn'  •  :,1  (  m‘ nnic  Usine  a  number  0!  l.ncari/ations  is\er  lite  lull  r.mgc  0.  s.d.e 

'ts  ,  ni  em.'i.t  a  niJ  Meld  more  a^^s.rate  s;f;v...ation  ressdts  ,\s  ’i.e  salse  'rans.tmiis 
t'  r  moe  m  pos.t. 'iw  dnlemnt  •more  l^.^-ur!:cr  dns  rr-’an  m  sons'. .nw  a  -n’d 
s  :‘-s  ;  ;  1'  :s  'dc'-n/ed  ti..;r  'i..'  annro.iJi  uni  inipr'  .e  tiie  per:  nin.nne  n  '>  ’d 


Plant  (  haracteri/atinn 

‘  :  'ns  1^.11"  i'S'  t  t.i.s  .i.;',,i.  ressiwn  1  s;;;e  '.sum  jinei;’':  its  usesl 

_ 1'..'  ,  .e  'll  s  rs.i.e  se'.er.i.  vdt'crettt  ire.p.ietissss  t^  snairutsri/s  '..c 

'■  '...s  ■  s  -I  ...  :t;e.. ..  I’s  ;  :  tins  rs  r\  moie  .M^ur.w'.  .nut  !..s  : 'i.:">  t  n.'s: 

;•  n  .  ■  .  os  ;-o  ns’m  s  nr  t  ■  In..',  s d i . w M: on-  me  1^1, in'  IX  nne  tiin 


u-;il  help  :o  ‘urther  ccnfirm  the  analytitai  theory  presented  and  more  accurately  hx  the 
cntpirically  obtaiited  constants.  1  his  work  can  be  done  with  a  spectrum  analyzer. 

'  Model  Supply  and  Return  Pressure  Dynamics 

In  the  existing  ntodel  supply  pressure  was  assumed  to  be  constant  and  return 
pressure  .v.is  assumed  to  be  zero.  Since  there  was  no  accumulator  in  the  supply  line  of 
tl'.e  sssieru,  under  't^Jy  the  supply  pressure  varied  slightly  with  demand.  Return 
pressure  aittMueli  Imv  was  not  zero.  Inclusion  in  the  model  of  these  elTects  may 
imrru  .  e  s.uuuation  accuracy, 
d  State  Space  Modeling 

In  order  to  oht.r.n  a  linear  model  of  the  system  under  study  some  simplifmg 
assun.piions  were  made.  One  of  these  assumptions  was  that  the  arbitrary  load  force 
on  the  piston  il  j  ■  was  zero.  This  may  have  been  an  oversimplification  and  may  be 
the  reason  'iie  digital  sirnulatien  does  not  accurately  predict  system  perform.ance. 
represents  iViction  Ior..es  uruch  are  non-linear  and  vary  with  piston  velocity  and 
siire^ti.'n,  .A  suite  space  representation,  where  Fj^  is  included  as  one  of  the  state 
\  ,i;'!a'’'.es,  nt  r.  ciimar.ate  this  prcbleni  and  improve  the  digital  model  performance. 

'  Influence  of  Temperature  on  Ssstem  Performance 

SMbiiry  dependence  on  temperature  was  observed  during  testing  of  the  digital 
s',  stem.  .Although  it  ic  known  that  higher  fluid  temperatures  decrease  viscous  damping 
imi  thus  hise  an  c'.crali  destabilizing  elTect  on  the  system,  other  less-obvious 
•empera'ure  eilccts  mav  aSo  be  important.  The  relative  insensitivity  to  temperiture 
'..nder  an  i!  -g  con'rol  is  an  important  issue  yet  to  be  resolved. 

0  Dat.i  Acquisition  Ssstem 

1  iie  present  IIP  microccmputer  system  has  the  capacity  to  be  used  as  a 
..  ntr  u.e-  .ind  d.ita  acouisuion  sxstem  .\11  hard’.varc  required  is  m  place,  but  time  did 
n  pe-i',,'  ■'nc  d.". el.,  pment  ol  soltwirc  to  implement  the  data  acquisition  capabilities 
:  Me  s',  uem  I  i'.i'  'x  ould  be  usel'ui  lor  further  studies  of  the  hxdraulic  pruer  element 
::.,i  :  r  1  i"  .snrk  perl,  rmed  t.'.  support  gradu.ite  classes  m  adsaiwed  controls. 


APPENDIX  A 

ANALOG  SIMULATION  CODE  (DSL) 


•k 


*  SIMULATION  OF  THE  VALVE  CONTROLLED  PISTON 

(VCP)  SERVOMECHANISM 

*  (LINEAR  ANALOG  CLOSED  LOOP  MODEL) 


*  (BLOCK  DIAGRAM  CODING) 


K  * 


EETAE 

BP 


CO? 

CIG 

CIP 

EMAK 

ERR 

ETAF 

FL 

G 

GA 

u 


IC 

IMAX 

IREQ 

KCE 

KCI 

51 


?LC 

PLl 

?LMAX 

PS 


=  ACTUATOR  EFFECTIVE  AREA,  IN**2 
=  EFFECTIVE  BULK  MODULUS,  PSI 
=  LOAD  DAMPING  CONSTANT,  LB-SEC/ IN. 

=  EFFECTIVE  CAPACITANCE  OF  TRAPPED  FLUID,  IN**3/PSI 
=  INPUT  GAIN,  VOLTS/ (INPUT  IN  IN) 

=  ACTU.ATOR  INTEPJ^AL  LEAKAGE  COEFFICIENT,  IN**3/SEC-PSI 
=  MAXIMUM  VOLTAGE  FOR  FEEDBACK  TRANSDUCER 
=  ERROR  SIGNAL  JUST  AFTER  SUMMER,  VOLTS 
=  ACTUATOR  FORCE  EFFICIENCY 
=  EXTERNAL  APPLIED  LOAD  FORCE,  LB 
=  LOAD  SPRING  CONSTANT,  LB/ IN. 

=  SERVOAHPLIFIER  OUTPUT  GAIN,  MA/VOLT 
=  FEEDBACK  GAIN,  VOLTS/ IN 

=  LIMITED  CURRENT  SUPPLIED  TO  SERVOVALVE,  MA 
=  LINEARIZATION  MIDPOINT  FOR  VALVE  CURRENT,  MA 
=  SERVO  VALVE  MAXIMUM  CURRENT,  MA 
=  LIMITED  CURRENT  REQUIRED,  MA 
=  EFFECTIVE  LEAKAGE  COEFFICIENT,  IN**3/SEC/FSI 
=  V.iiLVE  FLOW  PRESSURE  COEFFICIENT,  IN*’'3/SEC  PSI 
=  LINEARIZATION  CONSTANT  FOR  VALVE  FLOW  GAIN,  IN-^*3,SE( 
=  MODIFIED  VALVE  FLOW  GAIN,  IN^*3/SEC-MA-SQRT ( PSI ) 

=  TOTAL  ACCELARATED  MASS,  LBS*SEC**2/ IN 
=  LIMITED  LOAD  PRESSURE  DROP,  PSI 

=  LIMEARIZATICN  MIDPOINT  FOR  LOAD  PRESSURE  DROP,  PSI 
=  UNLIMITED  LOAD  PRESSURE  DROP,  PSI 

=  ;'.AXIMUM  PERMISSIBLE  LOAD  PRESSURE  DROP  FOR  MODEL,  PS: 
=  SUPPLY  PRESSURE,  PSI 
=  FLOW  FOLLOWING  SUMMER,  IN’^’"3/SEC 
=  LIMITED  VALVE  FLOW,  IN*’'3/SEC 


★ 


■k 

k 

k 

k 

k 


MA 


*  QLl  =  UNLIMITED  VALVE  FLOW,  IN**3/SEC 

QLHAX  =  MAXIMUM  PERMISSIBLE  VALVE  FLOW  FOR  MODEL,  IN**3/SEC 
VT  =  TOTAL  VOLUME  UNDER  PRESSURE,  IN«’^3 

*  X?  =  PISTCN/LOAD  POSITION,  IN 

*  XPD  =  PISTON/LOAD  SPEED,  IN/SEC 

*  XR  =  INPUT  OR  COILMANDED  PISTON/LOAD  POSITION,  IN 


■k 

k 

k 


k 


k 


kkkkkkkkTKkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

k 


k 

INITIAL 

I  =  0.0 
XP  =  0.0 
XPD  =  0.0 

k 


CONST 


k 


AP  =1.46, 

3P  =103.6, 
CIG  =  4.0, 

CIP  =  3.1E-03, 
FL  =0.0, 

G  =0.0, 

H  =  4,0, 


10 

0.0, 

KV 

= 

0.00352, 

HT 

= 

8.517E-02, 

PS 

= 

800.0, 

PLO 

= 

0.0, 

VT 

44.0, 

BETAE 

= 

3.0E04, 

GA  =  110.0, 
EMAX  =  12.0, 
ETAF  =1.0, 

IMAX  =  200.0, 
QLHAX  =  9.95606, 
PLHAX  =  400.0, 
KCI  =  6.0E-03 


PA.RAM  XR  =  1.7S5 


*  THIS  PART  OF  THE  CODE  DEFINES  SOME  OF  THE  CONSTANTS 

k 


KCE  =  KCI+CIP 

KQI  =  KV*(PS-PL0)*’^0.5 

CCP  =  VT/(4.C*BETAE) 

T1  =  CCP/KCE 
T2  =  (MT/BP) 


DERIVATIVE 


^  THIS  PART  OF  THE  CODE  IS  THE  PLANT  FROM  THE  ANALOG  BLOCK  DIAGRAM 


Q  =  I*KQI-XPD’'A? 

Q2=  Q/KCE 

FL=  REAL?L(0.0,T1 ,Q2) 

F  =  PL*ETAF'A?-FL 
rl=  F/EP 

XPD  =  REALPL(0.0,T2,F1) 

x?:=  intgrl;o.o,x?j) 


S 


XP  =LIMIT(-3.0,3.0,XP1) 


DYNAMIC 

•k 


^  THIS  FART  OF  THE  CODE  REPRESENTS  THE  SUMMER  AND  THE  AMPLIFIER  GAIN 


k 


ERR  =  ((XR*CIG)-(XF’*^H)) 

IF  (ERR. GT. 1.6. OR. ERR. LT. -1,8)  THEN 
I=2CQ.0*ABS(ERR)/ERR 

ELSE 


I=ERR*GA 


ENDIF 


CONTRL  FINTIM  =0.50 
SAVE  0.010,  XP 
FRI.NT  0.010,  XP 
END 

G.RAFH  (DS=TE,K613)  TIME (UN=  '  SEC ')  ,  ,XP(UN='IN') 
LABEL  LINEAR  CLOSED  LOOP  ANALOG  MODEL 


2 


LEGEND 

EXPERIMENTAL  DATA 
MODEL  DATA  ' 


0.00  0.05  0.10  0.15  0.::0  0.25  0.30  0.35  0.40  0.45  0.50 

TIME  (SECONDS) 

Figure  B.I  .Analog  System  Step  Response  (0.4088  in.). 


!  LEGEND 

!  =  EXPERIMENTAL  DATA 

"  MODEL  DATA 


0.00  0.05  0.10  0.15  0.20  0.25  0.30  0.35  0.40  0.45  0.50 

TIME  (SECONDS) 


Figure  B.2  .Analog  System  Step  Response  (0.635 


2'  ^ 


LEGEND 

c  EXPERIMENTAL  DATA 
o  MODEL  Data  “ 


0.00  O.Oo  O.IO  O.lo  0.20  0.25  0.00  0.G5  0.40  0.45  0.50 

TIME  (SECONDS) 

Figure  B.3  Analog  System  Step  Response  (0.8625  in.). 


LEGEND 

□  EXPERIMENTAL  DATA 
S  TVrODEL  data 


a.ooo.05  n.io  0.15  n.20  0.25  0.30  0.35  n.40  0.45  o.5o 

'TT'.ttr 

Figure  B.J  .Analog  Svstcm  Steo  Response  i  1.11  in.). 


LEGE\D 

□  EXPERIMENTAL  DATA 


U.UCO.Oo  0.10  0.15  O.QO  0.25  0.00  0.35  O.-tO  0.45  0.50 

TIME  (SECONDS) 

Figure  B.5  Analog  System  Step  Response  (1.32  in.). 


LEGEND 

□  EXPERIMENTAL  DATA 


0.000.05  0.10  0.15  0.20  0.25  0.30  0.35  0.40  0.45 

TIME  (SECONDS) 

Figure  B.6  Analog  System  Step  Response  (1.57  in.). 


APPENDIX  C 

MICROCOMPLTER  CONTROL  SYSTEM  ALGORITHM  (BASIC) 


10  REM  CLOSLOOP 
20  REM 

30  REM  THIS  PROGRAM  IS  FOR 
40  REM  CLOSED  LOOP  POSITION 
50  REM  CONTROL  OF  A 
60  REM  HYDRAULICALLY  ACTUATED 
70  REM  PISTON 
SO  REM  K=KPC 
90  REM 
100  K=.05 
lie  SETTIME  0,0 

120  ON  KEY#  1,"INTERUPT"  GOTO  130 
130  DISP  "WHAT  IS  DESIRED  POSITION  ?' 
140  DISP 

150  DISP  "INPUT  -3.0  TO  3.0  INCHES" 

160  INPUT  P 

170  CLEAR 

130  KEY  LABEL 

190  R=P*3.3333 

200  OUTPUT  723  ;"IP  8T" 

210  SEND  7  ;  UNL  MLA  T.4LK  23  SCG  1 
220  ENTER  7  ;  F 
230  E=(R-F)*K 
240  S=SGM(E) 

250  IF  .4BS(E)>10  THEN  E=S*10 
260  OUTPUT  723  ;  "OP  3"  ,E,"T" 

270  IF  ABS(E)>.125  THEN  GOTO  200 
230  OUTPUT  723  ;"IP  8T" 

290  SEND  7  ;  UNL  MLA  TALK  23  SCG  1 
300  ENTER  7  ;  F 

310  IF  .4BS(R-F;>.125  THEN  GOTO  230 
320  IF  TIMEM20  THEN  GOTO  340 
330  GOTO  280 
340  DISP 

350  CIS?  "PROGRAM  END" 

360  END 
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APPENDIX  D 

DIGITAL  SYSTEM  SIML  LATION  CODE  (DSL) 


3ETAE 

BP 

CC? 

^  c;g 

*  CIP 
^  EMAX 
ERR 

^  ETAF 


*  ?L 
^  PLO 
^  ?L1 

*  PLHAX 


SIHUIATIOM  OF  THE  VALVE  CONTROLLED  PISTON 
(VCP)  SERVOMECHANISM 
(LINEAR  DIGITAL  CLOSED  LOOP  MODEL) 


(BLOCK  DIAGRAM  CODING) 


*r7r******x**7t***7c**QL0SSARY*’^*’'**’'**’'**’'*’'****^’'*’'’'’'*"*’^***’ 

=  ACTUATOR  EFFECTIVE  AREA,  IM**2 
=  EFFECTIVE  BULK  MODULUS,  PSI 
=  LOAD  DAMPING  CONSTANT,  LB- SEC/ IN. 

=  EFFECTIVE  CAPACITANCE  OF  TRAPPED  FLUID,  IN**3/PSI 
=  INPUT  GAIN,  VOLTS/ (INPUT  IN  IN) 

=  ACTUATOR  INTERNAL  LEAKAGE  COEFFICIENT,  IN**3/SEC-PSI 
=  MAXIMUM  VOLTAGE  FOR  FEEDBACK  TRANSDUCER 
=  ERROR  SIGNAL  JUST  AFTER  SUMMER,  VOLTS 
=  ACTUATOR  FORCE  EFFICIENCY 
=  EXTERNAL  APPLIED  LOAD  FORCE,  LB 
=  LOAD  SPRING  CONSTANT,  LB/ IN. 

=  SERVOAMPLIFIER  OUTPUT  GAIN,  MA/VOLT 
=  FEEDBACK  GAIN,  VOLTS/IN 

=  LIMITED  CURRENT  SUPPLIED  TO  SERVOVALVE,  MA 
=  LINEARIZATION  MIDPOINT  FOR  VALVE  CURRENT,  MA 
=  SERVO  VALVE  MAXIMUM  CURRENT,  MA 
=  LIMITED  CURRENT  REQUIRED,  MA 
=  EFFECTIVE  LEAKAGE  COEFFICIENT,  IN**3/SEC/PSI 
=  VALVE  FLOW  PRESSURE  COEFFICIENT,  IN’'*3/SEC  PSI 
=  CONTROLLER  GAIN,  (DIMENSIONLESS) 

=  LIME.ARIZATION  CONSTANT  FOR  VALVE  FLOW  GAIN,  IN**3/SEC  HA 
=  MODIFIED  VALVE  FLOW  GAIN,  IN^*3/SEC-MA-SQRT(PSI ) 

=  TOTAL  ACCELARATED  MASS,  L3S*SEC**2/ IN 
=  LIMITED  LOAD  PRESSURE  DROP,  PSI 

=  LINEARIZATION  MIDPOINT  FOR  LOAD  PRESSURE  DROP,  PSI 
=  UNLIMITED  LOAD  PRESSURE  DROP,  PSI 

=  MAXIMUM  PERMISSIBLE  LOAD  PRESSURE  DROP  FOR  MODEL,  PSI 
=  SUPPLY  PRESSURE,  PSI 


FLOW  followi:;g  summer,  :n**3/sec 

LIMITED  VALVE  FLOW,  IM''''3/SEC 
•JMLIMITED  VALVE  FLOW,  ir;*’'3/SEC 

M.AMIMLM  PERMISSIBLE  VALVE  FLOW  FOR  MODEL,  IN*’^3/SEC 
TOTAL  VOLUME  UMDER  PRESSURE,  IN*"3 
?I5TCO;/LOAD  POSITION,  IN 
PISTON, 'LOAD  SPEED,  IN, /SEC 

INPUT  OR  COMMANDED  PISTON/LOAD  POSITION,  IN 


i^AAAAA-ArA'^ 


rA^ATvAArKKA 


INITIAL 


I  =  C.O 

:-:p  =  q  . : 

>:?D  =  0.0 


A? 

= 

1.46, 

10 

= 

0.0. 

GA 

= 

110.0  , 

3? 

= 

103.6, 

KV 

= 

0.C0352, 

EMAX 

= 

12. C  , 

DIG 

= 

3.3333, 

MT 

= 

S.5I7E-02, 

ETAF 

= 

1 .  c , 

CIP 

= 

S.1E-C3, 

PS 

= 

300.0, 

IMAM 

= 

200.0, 

r  • 

= 

0.0, 

PLO 

= 

O.G, 

QLMAX 

= 

9 .95606 

'j 

= 

C.O, 

VT 

= 

44.0, 

PLMAX 

= 

4C0.0, 

K?  C 

= 

3.3333, 

0.05 

BETAE 

3.0E04, 

KCI 

= 

6.  :'E-03 

O-I  >;r  =  3.Q0 

THIS  P.ART  OF  THE  CODE  DEFINES  SOME  OF  THE  CONSTANTS 


KCE  =  KCI+CIP 
KQI  =  KV’^(PS-?L0)’^*0.5 
CCP  =  VT/(4.0*BETAE) 
n  =  /"P/KCE 
T2  =  ..IT/BP) 


DERIVATIVE 

’■  THIS  PART  OF  THE  CODE  IS  THE  PLANT  FROM  THE  A.NALOG  BLOCK  DIAGRA 


Q  =  I*KQI-XPD*AP 
02=  Q/KCE 

PL=  REALPL(0.0,T1,02) 


F  =  PL-^ETAF^AP-FL 


Fl=  F/5P 

K?D  =  REALPL(0.0,T2,F1) 
;■:?:=  imtgrl(o.o,xpd) 

:•:?  =L:Mir(o.o,6.oo,xpi; 


THIS  PART  OF  THE  CODE  REPRESENTS  THE  MICROCOMPUTER  CONTROLLER  AND 

amplifier  gain 


ERR  =  ( fXR*CIG)-(XP*H))*KPC 
IF  (ERR.GT.1.3.0R.ERR.LT.-1.S)  THEN 
IREQ=200.0’^ABS(ERR)/ERR 

ELSE 


IREQ=ERR^GA 

END  IF 

Y=  DELAY(035Q,0.350,IREQ) 

X=  IMF'JLS(0. 0,0.350) 

I re; 2=  ZH0LD(X,Y) 

I  =  LIMIT (-IHAX,IMAX,IREQQ) 


CONTRL  FINTIM  =02.00,  DELT=0.00l 

IIEIHCD  RKSFK 

SAVE  0.010,  X?,I,ERR 

PRI.NT  O.CIC,  X?,I,ERR 

E!:d 

GRAPH  ;jE=TEX61S)  TIME (UN= ' SEC ') ,  XP(UN='IN'),  XPEX ( UN= ' IN ' ) 

LABEL  LIN'EAR  CLOSED  LOOP  DIGITAL  MODEL 

END 

STOP 
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APPENDIX  E 

DIGITAL  LOGIC  TEST  DATA 


RKSfX  INTEGRATION  METHOD  USED 

DSl  OUTPUT  LISTING,  GROUP  1 


TINE 

aP 

T 

ERR 

0  .ODIOOE+OO 

0  .COOOOE+00 

O.OOOOOE+OO 

0.49999 

1.000C0E'02 

O.OCOOOE+00 

O.OOOOOE+OO 

0.49999 

2 .00  :'00E-G2 

Q  ,OOOOOE+O0 

O.OOOOCE+00 

0.49999 

3 .0C000E-C2 

O.OOOOOE+00 

O.OOOOOE+CO 

0.49999 

4.00000E-02 

G.OOOOOE+OO 

0 .OCOOOEtOO 

0.49999 

5 .OCOOOE-02 

O.OOOGOE+OO 

O.OOOOOE+OO 

0.49999 

6.OC0:OE-O2 

O.OOOCOE-rOO 

O.OOOOOE+OO 

0.49999 

7.000CGE-C2 

0,OOOGOE+O0 

O.OOOOOE+OO 

0.49999 

3.00C0CE-G2 

G.OOOOOE+OO 

O.OOOOOE+OO 

0.49999 

9.CCC00E-02 

O.OCOOOE+00 

O.OOOOOE+OO 

0.49999 

1  .OOOOOE-C'l 

G.OOOOOE+OO 

O.OOOOOE+OO 

0.49999 

O.llCOO 

O.OOOOOE-OO 

O.OOOOCE+OO 

0.49999 

C.  12000 

O.OOOCOE-OO 

O.OOOOOE+OO 

0.49999 

0.13000 

0.00000E*00 

O.OOOOOE+OO 

0.49999 

0.14000 

O.OOOOOE+OO 

O.OOOOOE+OO 

0.49999 

0 . 1 D  0  0  0 

O.OOOCOE+00 

O.OOOOOE+OO 

0.49999 

O.loOOO 

O.OOOOOE+OO 

O.OOOOOE+OO 

0.49999 

0.17000 

O.OOOOOE+CO 

O.OOOOOE+OO 

0.49999 

0.13000 

O.OOOOOE+OO 

O.OOOOOE+OO 

0.49999 

0.19000 

O.OOOOOE+OO 

O.OOOOOE+OO 

0.49999 

0. 20  000 

O.OOOCOE+OO 

O.OOOCOE+OO 

0.49999 

0,21000 

O.COOOOE+OO 

O.OOOOOE+OO 

0.49999 

0.22000 

O.OOOOOE+OO 

O.OOOOOE+OO 

0.49999 

C  .23000 

O.OOOOOE+CO 

O.OOOOOE+OO 

0.49999 

0.24000 

O.OOOOOE+OO 

O.OOOOOE+OO 

0.49999 

0.25000 

O.OOOOOE+OO 

O.OOOOOE+OO 

0.49999 

0.23000 

O.OOOOOE+OO 

O.OOOOCE+OO 

0.49999 

C.27GG0 

O.OOOOOE+OO 

O.OOOOCE+OO 

0.49999 

0 . 23000 

O.OOOOOE+CO 

O.OOOOOE+OO 

0.49999 

C. 29000 

0 .  OCOGOE+OO 

O.COOOOE+OO 

0.49999 

0.30000 

O.OOOOOE+OO 

O.OOOOOE+OO 

0.49999 

0.31000 

O.COOOOE+QO 

O.OOOOOE+OO 

0.49999 

0 . 32CCC 

0 . OOOOOE+00 

O.OOOOOE+OO 

0.49999 

0.330C0 

O.OOOOOE+OO 

O.OOOOOE+OO 

0.49999 

0,34000 

O.OOOOOE+OO 

O.OOOOOE+OO 

0.49999 

0,35000 

O.OCOOOE-fOO 

54.999 

0.49999 

0 , 3&C00 

7 .01423E-03 

54.999 

0.49383 

3i,;o 

1,9727 

41 . 165 

C. 17121 

1.9594 

41.165 

0.16344 

z  3 1  r 

■>  '  ^  ^ 

41.165 

0.16E56 

j  -t  _ 

2.2227 

41.165 

C  .  1,6233 

3  5  3  > 

:.C393 

41 . 165 

0.16011 

~  C  V  - 

2 .2  560 

41.165 

C. 15733 

•  3 

2 .  :'27 

41.165 

0.15456 

44C0 
*1 ' 


J  J 

4  6  D  C 
47'::' 
4  c  I  '7 
4977 
::CG 
■  '  70 


•  «.  V  <J 
00 
:o 

5400 

:5:0 


0  JO 

5  3  0  C 


i : :  0 

jZOO 


^30  0 

•  0  ~  /y 


3  O’ 0  0 


c  r:  ,  j 

70  0  0 
7  1  '  0 


Mo; 


,0393 
.  1C  60 
.1226 


J75 


,  1498 
.  1612 
.1722 
.1833 
.1942 
.2051 
.2160 
.2263 
.2377 
.  2466 
.2595 
.2704 
2.2312 
2.2921 
2.3030 
.3139 
.3247 
.3356 
3465 
3574 
3633 
3791 


UO 


2.40C9 

2.4118 

2.4226 

2.4335 

2.4444 

2.4553 

2.4662 

2.4770 

2.4579 

2.4935 

2.5097 

2.5192 

2.52c3 


41.165 


41 . 165 
26.365 


.  385 


26.335 

26.885 

26.835 

26.885 

26.885 

26.835 

25.585 

26.385 

26.385 

26.535 

26.335 

26.385 

26.885 

26.665 

26.835 

26.385 

26.385 

26.885 

26.835 
26.335 

25 . 835 
26.385 
26 . 335 
26.885 
26.385 
26.385 
26.685 
26.385 

26 . 835 
26 . 385 
25.655 
26.335 
16 . 035 
16.035 
16.085 


.  1 Dx  .  O 

,  149CC 
,14623 
,14375 


0.14169 


0.12979 
0.13794 
0.13612 
C .12430 
0.13249 
0.13C67 
0.12386 
0.12705 
0.12522 
12342 
12161 
11979 
11798 
11617 
11435 
11254 
11073 
10892 


0.1G710 


0, 

0. 

0. 


10529 
10348 
10166 
9.93502E-02 
9 .30371E-02 
9 .62241E-Q2 
9.44110E-C2 
9.25979E-02 
07849E-02 
S9713E-02 
7158SE-02 
53457E-02 
8.35327E-02 
3. 17196E-C2 
8.01361E-02 
7 . 38649E-02 


tt 


.■’500 

2.5337 

16.055 

7.77112E-02 

■’90  0 

2 . 5404 

16.035 

7.66011E-Q2 

C  V  -  - 

2. 5470 

16.035 

7 . 55069E-02 

clOG 

2.5535 

16.035 

7.44137E-02 

s::o 

2.5600 

16 . 035 

7.33327E-02 

550  0 

2.5665 

16.085 

7.22475E-C2 

34-0 

2 . 5~3  0 

16.035 

7 . 11625E-C2 

5  ECO 

2.5795 

16.035 

7.00777E-C2 

2.5S60 

16.035 

6.39930E-02 

a'oc 

2.5925 

16.085 

6.79032E-C2 

35  00 

2.5991 

16.035 

6 . 66235E-02 

Z  Z, 

2.6056 

16.085 

6.57337E-a2 

500  0 

2.6121 

16.035 

6.4o540E-G2 

910  0 

2,6186 

16.065 

6 . 35693E-02 

2.6251 

16.035 

&.24845E-02 

9200 

2.6216 

16.085 

6.13993E-G2 

940  0 

2.6331 

16.065 

6.03150E-02 

9500 

2 . 6446 

16.085 

5.92203E-02 

96  00 

2.6511 

16 .065 

5.31456E-02 

9-7CG 

2.6576 

16 . 035 

5.70608E-02 

9  500 

2 .6641 

16.035 

5 . 59761E-02 

9  500 

2.6706 

16,035 

5.489:3E-C2 

0000 

2.6772 

16.035 

5.33C56E-02 

LIST  OF-  RLFERLNT'LS 


[  inch,  MB..  .\l:i:>oi'on;purer  CuniiiA  of  a  Uyciraulic  Power  Eienitni.  .M.L.  !  hesi<. 
Na-.al  Po'-tgraJuatc  School,  CA,  December  lhS6 

Merritt.  H  E:..  Ilyd>\iui:c  ('-■niroi  Systcois,  pp.  2-4.  84-81’.  141.  148.  .lohn  iie>  A 
Sons,  lire,, 

Henke,  R.  W ..  U  hy  doC'io  i  everyone  use  electronic  cmiirol.P .  Hydraulics  A 
Pneumatics,  pp,  51-54.  June  lhS5 

.Xstrom.  J.  K.  and  W’iitenmark  B.,  Computer  Controlled  Systems  I'/ia-ry  jrd 
Ceuyn.  p.  5,  Prentice-IIali  Inc..  I'^.84 

Rosko,  Joseph  S..  Diyitai  Simulation  of  Physical  Systems,  p.  106.  .AdJison-W'eslcy 
Pu'nlisinne  Co..  10“2 

Ogata.  Katsuiiiko  .'.["dern  Control  Enyineeriny,  pp.  625-6.' I.  Prentice-Hall  Inc., 
I'.ro 


IMTIAL  DISTRIBLTION  LIST 


Dc'.ense  Technical  Ini'orniation  Center 
La  moron  Saonon 
Alexandria.  \'A  22'('4-6145 

Lii'rai';.',  Code  iil42 
\a\al  Pcsteradiiatc  School 
\Ionterey,  CA 

Ciiairrmin,  Code  Iv 
Mecnanicai  Engineering  Department 
Na\ al  Postgraduate  School 
Mcntere}',  CA  43443 

Prolesscr  R.If.  Nunn.  Code  69\n 
Mecnanicai  Engineering  Department 
Na\al  Po'-tgraduate  Sclioci 
Moiierey.  C.\  43443 

Commandant  1  G-ENE.4fh 
L-S.  Coast  Guard  Meadquarters 
:ii  (I  2nd  and  V  st.  SAG 
M  ashington.  D.C.  2')i'')3 

c 0 rnma nda n t  ( G -  PT T  -  2 ; 

L.S.  Coa'it  Guard  Headquarters 
2  P'd  2nd  and  st.  S.W. 

A  ii'hington.  D.C.  doGG 

Gommandmg  OlTicer 
L  S.C.G.C,'.  -SPENCER  i AVMnC-905) 
c  0  G.G.  Support  Center 
42’  Comur.ercial  St. 

\L\  (J2109 

.Attn:  LI.  Ian  Grunther 


